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ABSTRACT 

We present a ready-to-use method to constrain the density distribution in early-type 
galaxy lenses. Assuming a power-law density profile, then joint use of the virial theorem 
and the lens equation yields simple formulae for the power-law index (or logarithmic 
density gradient). Any dependence on orbital anisotropy can be tightly constrained 
or even erased completely. Our results rely just on surface brightnesses and line-of- 
sight kinematics, making dcprojcction unnecessary. We revisit three systems that have 
already been examined in the literature (the Cosmic Horseshoe, the Jackpot and 
B1608+656) and provide our estimates. Finally, we show that the method yields a 
good approximation for the density profile even when the true profile is a broken 
power-law, albeit with a mild bias towards isothermality. 

Key words: Gravitational lensing: strong - Galaxies: elliptical and lenticular, cD - 
Galaxies: kinematics and dynamics - dark matter 
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1 INTRODUCTION 

When a galaxy acts as a gravitational lens, its mass can 
be probed independently via analyses of lensing properties 
and stellar kinematics. This combination has proved a pow- 
erful way of breaking both the mass-sheet degeneracy in 
lensing and the mass-anisotropy degeneracy in dynamics. 
It was first tackled by combining the lens equa tion with 
the Jeans equations (cf iTreu fe Koopmansl I2004T ) , usually 
with the assumption of a constant velocity anisotropy pro- 
file. Of course, Jeans analyses give no guarantee that the 
chosen models corre spond to physically consistent distri- 
butio n functions (cf lAn fc Evans! 120091 : ICiotti fc Morgantl 
2010). Latterly, sophisticated methods to build distribution 
functions of lensing galaxie s using the Schwarzschild met hod 
have been developed (e.g. iBarnabe fc Koopmansl [20071 ). In 
both cases, though, observational data are often not enough 
to constrain all of the parameters used. 

Here, we provide a simple method for combining lens- 
ing and dynamics. Our starting point is the virial theorem, 
which has the significant advantage that no assumption as to 
the anisotropy profile is ever needed. As a concrete case, we 
examine the scale-free density profiles p oc r~ 7 . Given an ac- 
curate photometric profile and an estimate of the Einstein 
radius from strong lensing, our method straightforwardly 
yields the value of the logarithmic gradient of the density 7. 
It does not require any photometric deprojection and it is 
robust under changes in orbital structure of the lens galaxy. 
By considering three benchmark lenses - the Cosmic Horse- 
shoe, the Jackpot and B1608+656 - we demonstrate that 
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our method performs very creditably against more sophisti- 
cated algorithms. Given its simplicity, we believe it provides 
the natural first step in combining lensing and dynamics. 



2 THE VIRIAL THEOREM AND THE LENS 
EQUATION 

One global constraint for galaxy lenses is given by the pro- 
jected mass within the Einstein radius. Another is given 
by the virial theorem, which links the average kinemat- 
ics of stars to the total mass density and luminosity. In 
other words, lensing probes projected masses, while pho- 
tometry and line-of-sight motions together trace the three- 
dimensional potential, so a combined analysis allows us to 
infer properties of the total mass profile. Here, we will quan- 
tify this framework by exploiting only the direct observables: 
surface brightnesses and projected kinematics. 

2.1 A Global Constraint 

The basic observables are the Einstein ring radius _Re, the 
surface brightness profile of the lens galaxy T,(R), together 
with the mean v\ os (R) and mean-square velocities w 2 i os (_R) 
of stars along the line of sight at projected radius R. 

Let us start with the lens equation. For a spherical lens, 
this gives the projected mass within the Einstein ring radius 
as 



M P (R E ) = 7rE cr 7? E 2 , 



(1) 



where S cr is the critical s urface density (see e.g., 
ISchneider. Ehlers fc Falcol 1 1992ft For a power-law lens with 
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an assumed total (i.e., luminous and dark matter) density 
p — po(r /ro) -7 , the projected mass within the Einstein ring 
M p (Re) reduces to a standard integral which can be done 
analytically. This yields: 



Po r o 



fl — K)E cr i?E 



7 -x (3- 7 )r( 7 /2) 
2V5FT(( T -l)/2) 



(2) 



where we have allowed for the presence of an ex ternal mean 
conve rgence k inside the Einstein radius (e.g., ISuvu et al I 
1201(f) . The external shear, even if important in modelling 
the lens, does not enter the mass normalization. 

The second ingredient is the virial theorem. For a spher- 
ical system, the virial theor em can be formulated in terms 
of the surface brightness as (|Agnello fc Evan"s"| [2012') 



(v 2 los)L = 



16ttG 
3Lt 



E(R)R 



p(r)r 2 



-.drdR (3) 



Here, Lt is the total luminosity, and we have introduced the 
notation 



{A)L = J- 

T Jo 



AT,(R)RdR. 



(4) 



to denote luminosity averages. 

For power-law density profiles, eqs ([2]) and Q give us 
immediately 
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(5) 



Here D s , D\ and D\ s are the angular-diameter distances to 
the source, to the lens and from lens to source respectively. 
The function T is analytic, namely 

r(x/2)r((5-z)/2) 



T(x) = 



(6) 



r(2-i/2)r((x-i)/2) • 

It is worth pausing to examine eq ([5j), which is the first 
main result of this Letter. Every quantity is known, or can 
be computed directly from the observations, with the sole 
exception of 7. Hence, the power-law index 7 depends just 
on global, measurable quantities. There is no dependence 
whatsoever on the orbital anisotropy. 

The magnitude of the exte rnal convergence k is affected 
by the mass-s heet degene r acy ( Falco et al.ll 19851 ) . Following 
the results of ISuvu et al I (|201C ~) . we adopt a prior distribu- 
tion for k with zero mean and ~ 0.05 dispersion. This trans- 
lates into a ss 5% uncertainty in the mass-normalization, 
or ~ 2.5% in the velocity dispersion. The latter is small 
compared to the typical w 15% statistical uncertainties of 
present-day measurements. 

2.2 Aperture Corrections 

In practice, the kinematics are measured inside a finite 
aperture radius i? ap and blurred by a point-spread func- 
tion (PSF). Hence, we re-state the results in terms of 
the aperture-averaged line of sight motions, introducing an 
aperture-correction: 



(v 2 los)l 



(7) 



(^ 2 io S )L«ii ap )(l + e) • 

The function e characterises how the aperture-average is re- 
lated to the luminosity-average over the whole system. If 
the kinematics are measured well beyond the effective ra- 
dius, or if the latter is small compared to the fiber aperture, 
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Figure 1. Logarithm of the aperture correction 1/(1 + e) as a 
function of 7. The lines are plotted for a de Vaucouleurs surface 
photometry and i? ap /i? e = 0.5,0.6,1,2,3,10. Smaller aperture 
radii correspond to steeper aperture corrections. 



then the correction is small and e as 0. For lenses in the 
SLACS samp le, R au is ~ 3/5 of the de Vaucouleurs effective 
radius R c (cf lAuger et all boich . Hence we need to inquire 
whether averaging within the aperture radius (rather than 
to infinity) introduces an effective dependence on orbital 
anisotropy. Let us assume, as a starting point, that the sys- 
tem is isotropic within the effective radius, as suggested by 
dynamical analyses of near by elliptical galaxies (jSaglia et al.l 
1 19931 : 1 Gerhard et al J [20011 ). We regard this as an informed 
prior on orbital structure and examine in a second step the 
consequences of anisotropy. We integrate the isotropic Jeans 
equations and reconstruct the kinematic profile to quantify 
how much of it is left out of the aperture-average. For a 5- 
function PSF, the result can be stated explicitly in terms of 
the surface brightness, namely 



1(0, 7) 



7(0,7) -I(R,p, 7) 



(8) 



were i<ij ap is the luminosity from projected radii less than 
the aperture radius R ap . We have also defined 



E(i?) J R 3 -^( 2/ /7?, 7 )di? 



(0) 



where the kernel g(y,j) is given by 



g{x,i) = (4-7) 
+ (7- l)x 2 



X 3—7 



da 



X 1 — 7 



du 



(10) 



Eq is the second main result of this Letter. Here, the 
aperture correction is expressed as a double quadrature that 
depends on the power-law density exponent 7. Fig. [T] shows 
the aperture correction 1/(1 + e) for a de Vaucouleurs pro- 
file with different values of R^p/Re- For lenses in the SLACS 
sample, 7 ~ 2, so that the correction is generally steeper in 7 
than the right-hand side of eq ([5]) and only weakly dependent 
on the details of photometry. This behaviour is understand- 
able, since the additional requirement of pressure isotropy 
reduces the freedom in 7 to yield physically consistent solu- 
tions. 

The aperture corrections also can be evaluated for 
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anisotropic velocity dispersions, such as 

/J(r) = 1-^=0,-^. (ii) 

The systematic error on 7 given by neglecting anisotropy is 
limited in modulus to < 5% for the SLACS sample, when 
P(R C ) varies between -0.5 and 0.5 and a de Vaucouleurs lu- 
minous profile is used. The correction to e is always negative 
for /3o > (radial anisotropy) and positive for /?o < (tan- 
ge ntial anisotropy). This i s an extension of the result found 
in iKoopmans et al.l (|2009l 'l. where Hernquist or Jafie lumi- 
nous densities with uniform anisotropy were used to quantify 
the effect. We conclude that at least as far as aperture cor- 
rections are concerned, velocity anisotropy is a secondary 
effect. 

2.3 Algorithm 

Eqs © and (JH]) provide a prediction for the root-mean- 
square velocities along the line of sight, for which it is con- 
venient to introduce the shorthand 

«( 7 ) = \A^los>L. (12) 

On the other hand, observations give us a value and uncer- 
tainties for the same observable, say a ± 5a. 

With larger fiber apertures, or extended long-slit mea- 
surements, there is no need to make any aperture correc- 
tions. So, for each value of 7, eq ((5]l yields the virial predic- 
tion and the likelihood on the density exponent is simply 

pMocexp-I^M-^. (13 ) 

If the kinematics are not measured over a sufficiently wide 
radial range, then both sides of eq (JSJ are divided by the 
aperture correction 1 + e, as defined in eq ((HJ, to get the new 
prediction for the line of sight velocity dispersions and the 
new likelihood for 7. This step is critical when the aperture 
is a small fraction of the effective radius, as 11(7) may be 
very flat or even non-monotonic (Fig. [2]); in this case, if the 
aperture correction is ignored then the inference on 7 may 
be biased and imprecise. 

As illustrated, the combination of lensing and dynam- 
ics has been reduced to the routine procedure of evaluat- 
ing a single non-linear function. When the PSF has a non- 
negligible width, the aperture-averages are more represen- 
tative of the luminosity averages. In that case, the aperture 
corrections are smaller than predicted by eq jSJ, which as- 
sumes a 5-function PSF. The density exponent must then lie 
between the uncorrected and aperture-corrected estimates. 



3 EXAMPLES 

We first discuss three particular lens systems, comparing 
the performance of our methocj]] with results in the litera- 
ture. The first quoted error bars from our method refer just 
to the statistical uncertainties propagated by line of sight 
kinematics. The sources of systematic error (anisotropy, ex- 
ternal convergence and PSF) are briefly touched upon case 

1 Code to implement our algorithm is available at 
http: / /www. ast . cam. ac .uk/~aagnello/virialLensing/ 
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Figure 2. The relationship between the model power-law index 7 
and the model velocity dispersion v without aperture corrections 
(dashed) and including aperture corrections (solid) for the 3 ex- 
ample lenses discussed in Section 3 (red, labelled H): Horseshoe; 
black, labelled J: Jackpot; blue, labelled B: B1608). In addition to 
large quantitative changes in the predicted dispersions, the aper- 
ture corrections also yield significant qualitative changes; in all 
cases the relationship gets steeper when the aperture correction 
is applied. This leads to more precise inference, particularly in 
the case of the Jackpot, where the flatness of the uncorrected re- 
lation implies little information on the slope given an observed 
dispersion, and the non-monotonic nature of ^(7) also leads to a 
(biased) bi-modal inference. 

by case and incorporated as a second uncertainty in the 
aperture-corrected estimates. We adopt as this systematic 
uncertainty 5% of the mean value, even though this might 
be too pessimistic. Finally, we examine the consequences of 
applying our single power-law inference to broken power-law 
density profiles. 

3.1 The Cosmic Horseshoe Lens (SDSS 
J1148+1930) 

The C osmic Horseshoe was discovered by iBelokurov et al.l 
|2007fl . It is a massive early-type lens at a redshift zi = 0.44 
surrounded by a loose group . The extended source galaxy at 
z s = 2.381 |Dve et al.ll2008T ) is lensed into a nearly complete 
Einstein ring. Early studies used estimates of the velocity 
dispersion deri ved from low re s olutio n spectra. 

Recently, ISpiniello et al.l l |201ll ) used high resolution 
spectra derived from X-Shooter on the Very Large Telescope 
to provide line of sight kinematic profiles out to the effective 
radius. Using the Jeans equations, they find a logarithmic 
density slope of 7 = 1.72±g |. We use their kinematic data 
(a = 352 ± 19 inside a circularised aperture of 0.84") with 
a de Vaucouleurs luminous profile (with an effective radius 
of 1.96"; the Einstein radius is 5.1"), to get an estimated 
value of 7 = 1.76 ± 0.08 without aperture corrections and 
1.80 ± 0.04 ± 0.09 when aperture co rrections are included . 
This is perfectly compatible with the lSpiniello et al.l (|201lh 
estimate within the uncertainties. It is worth remarking that 
their Jeans analysis is fit to the line-of-sight velocity disper- 
sion profile, neglecting any internal rotation. Our method 
exploits the whole velocity second moment {v 2 i ob )l, (i.e. 
pressure and rotation), which is slightly higher than the ve- 
locity dispersion alone. Moreover, we rely just on the surface 
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brightness profile E(i?) instead of using an approximation 
for the deprojected density, which they are forced to do so 
as to employ the Jeans equations. 

3.2 The Jackpot Lens (SDSSJ0946+1006) 



The Jackpot Lens was discove red by | Bolton et al.l (|2008T ) 
and subsequently analysed by iGavazzi et al. (|2008f ). The 
system is unusual in that it possesses two Einstein rings, 
as two different high redshift sources are lensed by the same 
foreground elliptical galaxy at a redshift of Z\ = 0.22. 

With the photometric, k inema tic, and lensi n g dat a pro- 
vided by lAueer et~ai1 l|2009l ) and lAueer et al.l ^PIOD . our 
aperture-corrected power-law exponent is 7 = 2.01 ±0.15 ± 
0. 1. This is perfectly consistent with the estimate 2.01±0.18 
of lAueer et al.l l)20ld ). where the uncerta i nties f rom the PSF 
are fully accounted for. ISonnenfeld et all (|2012h use a strong 
prior on 7 obtained by more detailed lensing analyses and in- 
formation from both sources to find 7 = 1.98±0.03, in good 
agreement with the previous estimates. Using their double 
Sersic photometry, our method yields 7 = 2.06 ± 0.04 ±0.1 
if we use the inner Einstein radius. On the other hand, 
their Jeans analysis approximates the luminous profile with 
two cored pseudo-Jaffe profiles, but it is not immediately 
clear how a flat-top luminous density can be supported 
self-consistently in a cusped potential like the scale-free 
ones , so the underlying dynamical model may be unphys- 
ical (|An fc EvansH2009T ). 

The value 7 = 2.20 ± 0.03 estimated by IVegetti et al] 
(|2010l ) is obtained by examining a small region around the 
inner Einst ein radius, whereas our method and the Jeans 
analysis in ISonnenfeld et al.l (|2012h provide an estimate of 
the global density exponent on the whole radial range out 
to the Einstein radii. This discrepancy hints that the true 
density profile is no t a simple power-law . One possibility is 
that the analysis bv lVegetti et al.l (|2010t ) probes the density 
profile where both baryons and dark matter are important, 
and the assumption of a single power-law is too facile. 

3.3 B1608+656 



B1608±656 (iMvers et~aT1 1 19951 ; ISnellen et al] 1 19951 ; 
iFassnacht et al. I Il996h consists of an early-type galaxy 
lens with a small companion at z\ = 0.63. The source is an 
active galaxy at z s = 1.394. 

In lKoopmans et al] |2003T ). the lensing morphology and 
time delays of the images were used to constrain the lensing 
potential. This was subs equently ref i ned b y incorporating 
kinematic information in ISuvu et al I (|2010l ). which allowed 
them to break the degeneracies between density slope and 
Hubble constant that would otherwise afflict the problem 
(| Wucknit j I2002T ) . Their complete analysis, assuming Jaffe 
or Hernquist lumino us densities, gives 7 = 2.08 ± 0.03. We 
use the same data as ISuvu et al I |2010;) (a = 260 ± 15 in a 
0.42" aperture, the effective radius is 0.6", and the Einstein 
radius of the main lens is 0.83") to obtain 7 = 2.39 ± 0.09 
for the uncorrected estimate and 7 = 2.14 ± 0.06 ± 0.11 
including aperture corrections. As with the Horseshoe, this 
estimate is consistent wi th but sligh t ly larg er than previous 
estimates, although the ISuvu et al I |2010l ) result is mostly 
driven by the lensing information from the extended source 
structure that is ignored in our analysis. 




Figure 3. The globally-estimated density exponent 7 g versus the 
inner power-law exponent 7 of the broken power-law model, as 
given in eq 1141 . The dashed band illustrates the model predic- 
tion when the aperture is very large (between 0.25 and 1 times 
the scale radius of the broken power-law) , while the grey band is 
for a small aperture that is 5 to 20 times smaller than the effec- 
tive radius. The widths of the bands demonstrate the effects of 
changing the size of the effective radius of the luminous tracer for 
fixed aperture and scale radius; for large effective radii the steep 
outer halo is given more weight and 7 g is more discrcpent from 
the true value 7. 



3.4 Extensions 

So far, we have assumed a power-law total density. Whether 
justifiable or not, this assumption is widely used out to 
length scales an order of magnitude gr eater than the effec- 
tive radius r B or Einstein radius -Re (cf lGavazzi et al]|2007l ; 
iHumphrev fc BuotefeOlOh . However, the normalisation from 
lensing is set by the projected mass within the (cylindri- 
cal) Einstein radius, which in turn includes mass at large 
(spherical polar) radii where the density profile may fall off 
differently from any inner power-law. 

As an example, let us consider the broken power-law 
density distribution 



p(r) = 



Po{r/r s y 



(l±r- 2 /rf)( 3 -x>/ 2 



(14) 



This behaves like p oc r ' for r < r s , while p — > r -3 for 
r>r s . The appearance of an explicit scale-radius r s means 
that the lensing and virial equations are not enough to fix 
all the independent parameters. The virial constraint is then 
just another piece of information to be incorporated into the 
lensing decomposition, albeit now with more complicated 
kernels than in eqs ([5} and ([8]). 

However, we may ask if approximating this profile with 
a global power-law p = pi(r/ri)~ ls gives a reasonable rep- 
resentation of the inner parts of the system. The virial esti- 
mates can help us answer this question. In particular, we find 
the global power-law index 7 g by fixing Re and (u 2 i os ) l<r^ 
to the ones given by the double-power law profile of eq (I14[) . 

Figure [3] shows bands of solutions corresponding to 
the physically plausible range of halo scale lengths with 
5 ^ r s /R-E < 20. The dashed curve represents situations 
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when the aperture is large (here, 5 times the Einstein ra- 
dius, or 0.25 to 1 times the scale radius). It is clear that for 
such large apertures the power-law approximation for the 
density profile is inadequate, with the central slope being 
grossly over-estimated for small values of 7 and systemati- 
cally under-estimated for large values of gamma. This trend 
persists when the aperture is small (the grey curve, where 
the aperture is equal to the effective radius, i.e., 5 to 20 
times smaller than the scale radius) but with a significantly 
smaller bias. 

The bias of the estimator 7 g is always towards isother- 
mality and in the case of small apertures (e.g., consis- 
tent with the apertures used for observations) the bias 
tends to favour 7 g ~ 2.1. This is remarkably close to 
the mean power-law in d ex inferred for the SLACS lenses 
ilKoopmans et al.l 120091 : lAuger et all |2010| ; iBarnabe et~al] 
l201ll ), although power-law indices inf erred without the us e 
of dynamics yield similar results fe.g.. |Pobke fc King||2006T ). 
Nevertheless, the magnitude of the bias is also smallest at 
7g ~ 2.1 and suggests that the power-law index estimates 
from lensing and dynamics robustly describe the central den- 
sity profile even if the true mass density distribution is a 
broken power-law. 



4 CONCLUSIONS 

We have presented a simple method to combine strong lens- 
ing and stellar dynamics without any need of prescribing the 
poorly known velocity dispersion anisotropy. Using measure- 
ments of the Einstein ring radius, the photometric profile 
and the global kinematics, we have shown how to obtain 
the logarithmic slope of the density 7. When the kinematics 
are measured on a small aperture, then the density slope 
7 and the velocity anisotropy /3 are intertwined, but proce- 
dures have been developed to characterise the influence of 
anisotropy robustly. Above all, our algorithm has the sig- 
nificant advantage of simplicity - all that is needed is the 
numerical solution of a single non-linear equation to obtain 
an estimate for the density slope 7. We have verified that 
our algorithm works well by comparison with analyses in the 
literature for three well-studied lenses. 

The set of hypotheses used here is common to a ma- 
jority of Jeans equation analyses, but the lack of any need 
for marginalisation over velocity anisotropy j3 in our method 
allows us to separate the statistical uncertainties from sys- 
tematics. Also, our formulae require only a profile for the 
surface brightness, avoiding any additional uncertainty from 
deprojection inherent in applications of the Jeans equations. 

Why does our simple method work so well? A major 
reason is that velocity anisotropy is never significant at the 
centres of early-type galaxies. If the kinematics are known 
within a finite aperture radius, then anisotropy can affect 
our result. However, this correction is always small, and 
isotropic models give a good guide as to its magnitude. The 
fact that the virial theorem is independent of the velocity 
anisotropy enables us to sidestep difficulties that Jeans or 
phase-space modelling must necessarily confront. 

The main limitation that this method shares with Jeans 
analysis is observational: aperture-averaged line of sight mo- 
tions must be suitably dealt with to account for different 
possible sources of bias. However, when care is taken in as- 



sessing these effects, our procedure gives the density slope 
7 in terms of directly observable quantities and can be used 
as a robust and quick method of carrying out lensing plus 
dynamics decompositions. 



AA thanks the Science and Technology Facility Council and 
the Isaac Newton Trust for the award of a studentship. We 
wish to thank Nicola C. Amorisco, Iulia T. Simion and Vasily 
Belokurov for significant feedback on the manuscript, as well 
as Alessandro Sonnenfeld for discussing some details of his 
analysis with us. The anonymous referee gave useful com- 
ments that helped us improve the paper. 



REFERENCES 

Agnello, A. & Evans, N.W. 2012, ApJL, 754, L39 
An, J. H., & Evans, N.W. 2009, ApJ, 701, 1500 
Auger, M. W., Treu, T., Bolton, A. S., et al. 2009, ApJ, 
705, 1099 

Auger, M. W., Treu, T., Bolton, A. S., et al. 2010, ApJ, 
724, 511 

Barnabe, M., & Koopmans, L. V. E. 2007, ApJ, 666, 726 
Barnabe, M., Czoske, O., Koopmans, L. V. E., Treu, T., & 

Bolton, A. S. 2011, MNRAS, 415, 2215 
Belokurov, V., Evans, N. W., Moiseev, A., et al. 2007, 

ApJL, 671, L9 

Bolton, A. S., Buries, S., Koopmans, L. V. E., et al. 2008, 

ApJ, 682, 964 
Ciotti, L., & Morganti, L. 2010, MNRAS, 408, 1070 
Dobke, B. M., & King, L. J. 2006, AA, 460, 647 
Dye, S., Evans, N. W., Belokurov, V., Warren, S. J., & 

Hewett, P. 2008, MNRAS, 388, 384 
Falco, E. E., Gorenstein, M. V., & Shapiro, 1. 1. 1985, ApJL, 

289, LI 

Fassnacht, C. D., Womble, D. S., Neugebauer, G., et al. 

1996, ApJL, 460, L103 
Gavazzi, R., Treu, T., Rhodes, J. D., et al. 2007, ApJ, 667, 

176 

Gavazzi, R., Treu, T., Koopmans, L. V. E., et al. 2008, 

ApJ, 677, 1046 
Gerhard, O., Kronawitter, A., Saglia, R. P., & Bender, R. 

2001, AJ, 121, 1936 
Humphrey, P. J., & Buote, D. A. 2010, MNRAS, 403, 2143 
Koopmans, L. V. E., Treu, T., Fassnacht, C. D., Blandford, 

R. D., & Surpi, G. 2003, ApJ, 599, 70 
Koopmans, L. V. E., Bolton, A., Treu, T., et al. 2009, 

ApJL, 703, L51 
Myers, S. T., Fassnacht, C. D., Djorgovski, S. G., et al. 

1995, ApJL, 447, L5 
Saglia, R. P., Bertin, G., Bertola, F., et al. 1993, ApJ, 403, 

567 

Schneider P., Ehlers J., Falco E.E., 1992, Gravitational 

Lenses, Springer Verlag, New York 
Snellen, I. A. G., de Bruyn, A. G., Schilizzi, R. T., Miley, 

G. K., & Myers, S. T. 1995, ApJL, 447, L9 
Sonnenfeld, A., Treu, T., Gavazzi, R., et al. 2012, ApJ, 752, 

163 

Spiniello, C, Koopmans, L. V. E., Trager, S. C, Czoske, 
O., & Treu, T. 2011, MNRAS, 417, 3000 



© 0000 RAS, MNRAS 000, [Lj{6] 



6 Agnello, Auger & Evans 

Suyu, S. H., et al. 2010 ApJ, 711, 201 

Treu, T., & Koopmans, L. V. E. 2004, ApJ, 611, 739 

Vegetti, S., Koopmans, L. V. E., Bolton, A., Treu, T., & 

Gavazzi, R. 2010, MNRAS, 408, 1969 
Wucknitz, O. 2002, MNRAS, 332, 951 



© 0000 RAS, MNRAS 000, (TJU 



